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Abstract 

In this paper we develop an approximate 
computational framework for simulation of the 
fluctuating flowfield associated with the complex vortex 
system seen at the side edge of a flap in a multi-element 
high-lift airfoil system. The eventual goal of these 
simulations is to provide an estimate of the spectral 
content of these fluctuations, in order that the spectrum 
of the noise generated by such flowfields may be 
estimated. Results from simulations utilizing this 
computational framework are shown. 

Introduction 

Sound from an aircraft induced purely by 
airflow not related to the engine is known as airframe 
noise. During approach its levels rival that of the 
engine, causing a threat to the successful certification of 
future subsonic aircraft. 

NASA’s Noise Reduction Program began a 
recent effort to study Airframe Noise in 1995, 
partnering with United States major airframe industries. 
NASA Langley’s role is to determine fundamental noise 
source mechanisms by relating sound generation to 
fundamental fluid mechanics. It is important to realize 
that airframe noise prediction methods currently 
employed by industrial designers are based partly on 
broad-brush scaling estimates and partly on empirical 
data; little if any direct information regarding the actual 
noise generation mechanism is used, and for the most 
part, the details of these generation mechanisms remain 
unknown. Our program involves building block 
experiments coupled to large scale and flight tests 
performed at NASA Ames and industry. The work to 
be described below details one aspect of this effort 
pertaining to advanced computational tools, in 
particular, the tools utilized to understand the unsteady 
fluid dynamics associated with the high-lift system. 


Prediction of Flow-Induced Noise 

Although the subject of considerable study for 
the case of jet flows [1], the direct prediction of the 
noise produced by locally-separated flows on aircraft 
components through the use of numerical simulation has 
received comparatively little attention. Additionally, 
such flowfields tend to be rather complex in their mean, 
thus experimental investigations of the fluctuations 
occurring in these flows are also lacking. Both of these 
shortfalls result from the spatial and temporal resolution 
that would be required for a complete description of the 
mean and fluctuating flowfield. For example, in a 
companion study to this article [2], steady Reynolds- 
averaged Navier-Stokes (RANS) calculations are 
compared with experimental results for a simple part- 
span flap model of a high-lift system. The computation, 
using a well-established code, required a mesh of over 
19 million grid points for adequate spatial resolution of 
the flowfield; these results are used in this study as a 
mean state for the computation of an approximate 
fluctuating flowfield, as will be described. Although 
this large computation, requiring about 50 hours of Cray 
C-90 CPU time, provides good spatial resolution of the 
steady flowfield, it is completely inadequate to 
accurately capture the majority of the fluctuations that 
are believed to generate noise in this flow. 

Additionally, these fluctuations are known to be 
broadband in frequency. Since the computational effort 
required to simulate unsteady phenomena is roughly 
proportional to the ratio of the highest- to- lowest 
frequencies, it is clear that the use of unsteady RANS to 
simulate these fluctuations is out of the question, at least 
for the near future, even if issues such as the calibration 
of the turbulence model for capturing unsteady flows 
could be answered. 

The experimental investigation with which 
these steady RANS computations are compared in [2] 
required no less effort to achieve an adequate 
description of the steady part-span flap flowfield [3]. 


1 This paper is declared a work of the U.S. Government 
and is not subject to copyright protection in the United States. 
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Two techniques primarily provided this data: pressure- 
sensitive paint and five-hole probe surveys. Data from 
the latter will be shown in a later section. 

Thus, it is clear that in order to provide some 
detailed estimate of the origin and frequency content of 
fluctuations in a complex aerodynamic flowfield such as 
that occurring about the high-lift system of a 
commercial subsonic transport aircraft, let alone 
directly predicting the noise generated by such a 
flowfield, some approximation is required. In this 
study, three major approximations are used to reduce 
the problem to a level that allows at least a rough 
estimate of the frequency content and directivity of the 
noise generated by the flow at the side edge of the flap 
of a high-lift system. Although approximate, these 
estimates are based on simulations of the true physical 
phenomena which generate the noise, and are of 
reasonable computational cost that parameter studies 
are possible. 

This first approximation comes in the 
invocation of the Lighthill acoustic analogy [4]. As will 
be seen in the following section, this allows the 
combined noise generation and propagation problem, as 
would be solved in a so-called “direct computational 
aero-acoustics” simulation using the unsteady 
compressible Navier-Stokes equations, to be divided 
into a computation of the fluctuations in the near field 
and a separate computation of the generation and 
propagation of the noise. The former computation is 
still somewhat too large to be manageable, due to the 
resolution requirements discussed above, requiring 
further approximation. The second approximation lies 
in the use of the incompressible Navier-Stokes 
equations to simulate the near-field fluctuations. As the 
local Mach numbers seen in the flap-edge flowfield are 
generally below 0.3 for typical applications, this is 
considered reasonable, but the reduction in 
computational effort afforded by this approximation is 
still less than required to allow for a full simulation of 
the near-field fluctuations in this complex flow. In this 
study, therefore, the use of the “temporal” 
approximation is used to reduce the simulation of the 
full three-dimensional flowfield into a series of two- 
dimensional simulations; this approximation will be 
discussed and justified in later sections. 

The use of local, approximate-but-rational 
numerical simulations to estimate the fluctuations in a 
flowfield that are expected to generate noise represents 
an improvement over previous studies, in that the 
physics of the generation of the fluctuations is not 
assumed a priori (beyond assuming that it is 
represented by the Navier-Stokes equations). In fact, 
the overarching philosophy embodied in the program of 
which this study is a small part (the NASA Airframe 
Noise Reduction subelement of the Advanced Subsonic 


Transport Program) can be employed in the prediction 
of any flow-induced noise. First, the steady flowfield 
must be understood, through a combination of detailed 
experiment and careful application of configuration 
RANS computations, using somewhat simplified 
building-block geometries. Next, the flowfield is 
examined for features that are capable of producing 
large-amplitude organized fluctuations of the proper 
scale and frequency, and approximate simulations to 
describe these fluctuations are carried out. Such 
fluctuations would reasonably be expected to be the 
primary noise generators in a flow, and would likely 
(but not necessarily) result from fundamental 
instabilities of the flowfield. It is well recognized that 
inflectional instabilities can occur even in nominally 
turbulent flows, such as shear layers, and result in the 
reorganization of steady mean vorticity into fluctuating 
vorticity of possibly large amplitude. The nonlinear 
interaction of these fluctuations, and/or the evolution of 
these fluctuations in a rapidly-varying mean flow, 
would be the mechanisms for noise generation. 


LighthilPs Acoustic Analogy 

What follows is a very brief and incomplete 
description of the fundamentals of Lighthill’s acoustic 
analogy, in order to illustrate how the theory allows, at 
least conceptually, for the noise source and propagation 
mechanisms to be separated. The reader is referred to 
the original reference [4] and works by Ffowcs- 
Williams [5] and Lilley [6] for more complete 
descriptions. 

The basis of the theory is an exact expression 
of the compressible Navier-Stokes equations. In the 
absence of boundaries and mean flow, this equation is: 



The LHS of the equation is clearly a linear wave 
equation, and depicts the propagation of density 
fluctuations due to sound. The RHS, which results from 
the nonlinear convection terms in the momentum 
equation, acts as a source to the medium outside the 
region of fluctuating flow. This source is assumed 
known, derived from solutions or estimates of the 
nearfield fluctuating flow. 

In the presence of (stationary) boundaries, an 
additional, surface source term appears [5]: 


dx, 


Pi m 


df_ 

dx ; 
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where f=0 defines the surface, and p i} contains the 
fluctuating surface pressure and viscous stresses. By 
asymptotic analysis, it may be shown that for smooth 
surfaces, this source term dominates. However, in the 
presence of comer or edges, current investigations in 
which full CAA computations are compared with 
acoustic analogy results indicate that the volume term 
must be included. 

Flap-Edge Flowfield 

Based on both experimental and computational 
studies [2] [3], a fair understanding of the steady 
flowfield in the vicinity of the edge of a flap of a high- 
lift system has been gained. Clearly, the difference in 
lift between flapped and unflapped sections of a wing 
with a part-span flap will result in a trailing vortex 
emanating from near the flap edge, as shown in Fig 1 
which depict the streamlines wrapping around the flap 
side edge. The details of the development of this vortex 
are surprisingly complex. Shown in Fig. 2 are contours 
of a quantity that approximates the streamwise 
component of vorticity, displayed in planes normal to 
the flap edge and chord line [3], denoted “crossflow” 
planes. As can be seen, the pressure difference between 
the upper and lower surfaces of the flap creates flow 



Figure 1 Streamlines in flap-edge flowfield 

Numerical Simulation Algorithm 

Assuming that it is possible to justify 
performing simulations in the crossflow plane only (as 
will be discussed in the following section), we wish to 
compute unsteady incompressible Navier-Stokes 
solutions in a geometry that includes the rectilinear end 
of the flap. A robust and highly-accurate algorithm for 
such solutions is described in [7]; the physical domain 
is divided into rectangular sub-domains as required, and 


around the edge. Two separation bubbles, with 
associated streamwise vorticity and rollup, are created 
at the upper and lower comers of the flap edge. The 
reattachment point of the side-edge vortex moves up the 
edge as the flow progresses down the flap, eventually 
reaching the upper comer. The side-edge vortex then 
travels over the upper comer, interacting and eventually 
merging with the upper-surface vortex. This leaves a 
single trailing vortex, which is continually fed with 
vorticity from the cylindrical shear layer that emanates 
from the lower edge comer. This mechanism of 
continual feed of vorticity into the vortex produces a 
strong jet-like flow in the core of the vortex, where 
streamwise velocities of over twice the freestream 
velocity have been measured. For this flap section, 
which was chosen for building-block study due to its 
simple geometry, the trailing vortex leaves the flap 
surface and is more than a vortex-diameter clear of the 
surface when it reaches the trailing-edge location. For 
other, more representative flap sections, careful RANS 
computations indicate that the vortex may remain in the 
vicinity of the surface or upper comer until reaching the 
trailing edge. This has implication in the amount of 
noise generated by fluctuations that develop in the 
vortex/shear-layer system. 



Figure 2 Contours of streamwise vorticity at flap 
edge 

tensor-product Chebyshev spectral collocation of the 
pertinent equations is employed in each sub-domain. In 
Fig. 3 is shown an example of such a multi-domain 
discretization. Note that the comers of the geometry are 
isolated at comers of sub-domains. Thus, the geometric 
singularities associated with the comers do not enter 
into the discretization, and do not disrupt the expected 
spectral accuracy. Two additional sets of sub-domain 
divisions are employed in the discretization in Fig. 3, 
which set up layers of sub-domains both above and to 
the right of the main computational region of interest. 
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In these outer sub-domains, the extensively-used buffer 
domain technique [8] is employed to provide a non- 
reflecting numerical outflow boundary condition which 
allows passage of fluctuations that reach the boundary 
of the overall computational region, preventing 
contamination of the interior solution. 

The “interface” conditions, the discrete 
equations which are used to couple the solutions across 
coinciding sub-domain boundaries, are specific to the 
particular equations being solved; for second-order 
equations, strong enforcement of C 1 continuity has been 
found to perform well, as well as various flux-balance 
conditions [7]. 



Figure 3 Multi-domain discretization of crossflow 
plane at flap edge 


A simple time-splitting algorithm is used to 
advance the incompressible Navier-Stokes equations in 
time. This robust and accurate method, described in 
detail in [9], is easily employed in this multi-domain 
context, as it reduces the solution of the Navier-Stokes 
equations to solution of a sequence of Helmholtz and 
Poisson equations. Solutions of these equations fit 
easily in the above-mentioned interface scheme, and are 
carried out very rapidly on a computer using a variant 
of tensor-product diagonalization [8], The resulting 
simulation code is so efficient that all of the results to 
be shown in this paper were obtained on an SGI Indigo 2 
desktop workstation, and results have been obtained in 
the past using a 133-MHz Pentium PC. 


Temporal Simulation Technique 

Consider a fluctuating disturbance with 
streamwise scale O(e), evolving in a mean flowfield 
which changes in the streamwise direction on a 0(1) 
length scale. This disturbance additionally has an 
0(1 /e) frequency, and convects at some 0(1) fraction of 
the mean velocity. For this situation, it is clear that a 
multiple-scales analysis is justified; the disturbance 
undergoes 0(1) changes in the time and distance for 
which the mean state changes by only 0(e). This 
formalism has been extensively employed in the area of 
simulations of transition to turbulence in boundary 
layers. 

Another way to look at the temporal simulation 
framework is to consider a computational frame that 
travels downstream at a constant speed. In this frame, 
the mean state appears to change on a long, 0(1) time 
scale, whereas the disturbance, due to its 0(1 /e) 
frequency, evolves on an O(e) time scale. In this case, 
the mean state may be considered as frozen. The 
implication in regard to the physics of the disturbance is 
that gradients that matter to the development of the 
disturbance lie only in the transverse, and not 
streamwise, direction. 

To demonstrate that the above ansatz of 
interchanging the roles of streamwise distance and time 
is justifiable in this flow, we performed a simulation 
which started with the potential flow solution around the 
edge and, while maintaining the potential-flow values 
on all outer boundaries, enforced no-slip on the surface 
of the flap. A modest Reynolds number of 1000 based 
on flap thickness was used. As seen in Fig. 4, the 
flowfield evolves in time with features strikingly similar 
to those seen in the actual flap-edge flowfield, looking 
in planes going from the leading edge to trailing edge. 
For instance, in frame b the development of the 
separation bubbles from the upper and lower comers 
may be seen, followed by the reattachment of the side- 
edge bubble traveling up the side in c, the side-edge 
vortex moving over the upper comer and merging with 
the top vortex in d and e, and finally the development of 
a single vortex fed by a cylindrical shear layer 
emanating from the lower comer in f. This sequence 
may be compared with the individual crossflow planes 
shown in Fig. 2. In fact comparisons (not shown) 
between the final frame and a crossflow plane taken 
from the RANS computation of Fig. 1 at approximately 
50% chord show surprising resemblance in overall 
structure and detail. 
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Figure 4 Streamlines of temporally-evolving flowfield 


The above indicates that the concept of 
equivalence of streamwise distance and time is at least 
heuristically justified, and within that framework, a time 
scale exists for the evolution of what may be considered 
as the mean state. A farther justification of the 
temporal framework, that disturbances occur within this 


mean state and evolve on a time scale much shorter than 
that of the mean state evolution, is observed when the 
Reynolds number in the above simulation is raised. As 
seen in Fig. 5, which shows snapshots taken at closely- 
spaced time intervals, small-scale vortical disturbances 
spontaneously develop and convect rapidly within the 
slowly-evolving overall flowfield. This is strong 
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indication that the vortex / cylindrical shear layer start was required for its initiation, and that the time 

system is indeed strongly unstable, since no “seeding” scale of these dominant instabilities is small compared 

of the disturbance but that remaining from the impulsive with that of the evolution of the mean flow. 



Figure 5 Spontaneous development of disturbances in temporally-developing flowfield 


From this we conclude that, to some level of 
approximation, the evolution of the disturbances 
resulting from the fundamental instability of the vortex / 
cylindrical shear layer system may be considered within 
a frozen mean state. This is a significant step toward 
the analysis of flowfields that are truly representative of 
those near the edge of flaps in high-lift systems; we may 
extract crossflow “cuts” from high-resolution steady 
RANS computations, and perform simulations using 
these solutions as the mean flowfield within which the 
disturbances that are conjectured to create noise will 
evolve. Results from such analysis are shown in the 
following section. 


Results 

In [2] is described the results from a series of 
high-resolution steady RANS computations, for flow 
about a generic part-span flap configuration which has 
been extensively used in a number of experimental 
investigations of flap-edge noise [10][1 1][3]. We 
utilize those solutions here as mean flowfields for 
disturbance simulations, interpolating the RANS 
solutions onto the spectral multi -domain discretization 
topology described earlier. These interpolations are 
carried out in planes that are taken to be roughly normal 
to the axis of the vortex system. In order to study the 
growth of disturbances of particular frequencies, a small 
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amount of local forcing is applied in these simulations, 
at a point on the lower surface of the flap that is a 
significant distance from the lower comer. The forcing 
is applied as a short region of oscillatory suction and 
blowing, and careful examination of the unsteady 
flowfield showed little if any influence of location and 
size of the forcing region on the overall results. This is 
an indication of the strength of the instability 
mechanism at work in this flowfield, and leads us to 
conclude that the details of what provides the initial 
energy to these fluctuations in the actual flowfield is 
probably irrelevant; what matters is the amount of 
amplification that the inflectionally-unstable flowfield 
provides, as a function of frequency. The flowfield 
details that influence this instability, the strength, 
location, and thickness of the cylindrical shear layer, for 
instance, are functions of the configuration and loading. 


Comparisons of simulation results with surface 
fluctuating pressure and farfield noise data for the 
experiment of [1 1] are forthcoming. 

For a plane located at a streamwise station of 
roughly 10%-flap chord, the mean flowfield consists of 
distinct side- and top-comer vortices, with the side-edge 
reattachment point at about mid-thickness. The shear 
layers over these vortices are relatively thin, so higher 
frequency disturbances are expected to dominate. In 
Fig. 6 is shown contours of the mean vorticity, and 
snapshots of the disturbance vorticity for forcing at 
three different frequencies. These frequencies are taken 
in the model scale, and the contours displayed are all 
normalized with respect to the maximum from each 
simulation. Each simulation required about 4 hours of 
SGI workstation time to reach statistical steady state. 



Figure 6 Contours of mean and disturbance vorticity, 10%-flap chord station 


Near the lower comer, fluctuation strength as 
measured by disturbance vorticity magnitude is greatest 
for the 15 kHz forcing, with 30 and 5 kHz disturbances 
being 30 and 39% smaller, respectively. At a point in 
the shear layer near mid-thickness, the lowest (5 kHz) 


disturbance dominates, with 15 and 30 kHz disturbances 
59 and 62% smaller. 

In a plane at approximately 50% flap chord, 
the side-edge and top vortices have merged, and the 
cylindrical shear layer / vortex system is well 
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established. Contours of mean and disturbance vorticity 
are shown in Fig. 7, again with three model frequencies 
of 5, 15, and 30 kHz forcing. The instability of the 
cylindrical shear layer is quite apparent, with the 5 kHz 
disturbances persisting with significant magnitude even 


as they are convected over the vortex. Maximum 
disturbance vorticity magnitudes in the shear layer are 
roughly equal for the three frequencies, demonstrating 
the broadband nature of the instability. 






Figure 7 Contours of mean and disturbance vorticity, 50% -flap chord station 


As mentioned above in the discussion of the 
Lighthill acoustic analogy, the fluctuating surface 
pressure is of interest in regard to the prediction of noise 
produced by a fluctuating flowfield. Contours of 
logioflP’suJ) are shown in Fig. 8, displayed with the time 
coordinate as the third spatial dimension. Apparent in 
the figure is the change in location of concentration of 


fluctuating surface pressure moving from side edge to 
top surface with frequency decreasing; this is expected 
from the disturbance vorticity contours of Fig. 7. It is 
believed that concentration of fluctuating surface 
pressure near geometric singularities such as comers is 
of particular importance in noise generation. 
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Figure 8 Contours of log^flp’.J) 


The use of broadband forcing in these simulations gives 
an idea as to the overall spectrum of disturbance growth 
in these flowfields. In Fig. 9 are shown contours of 
mean and disturbance vorticity for this same 50%-chord 
plane, in which broadband forcing is used. The forcing 
is constructed of 500 discrete frequencies between 700 


Hz and 48 kHz, with random initial phases. To measure 
the response spectrum of the flow, numerical probes are 
placed in the flow and on the surface, at which time 
histories of disturbance vorticity and fluctuating surface 
pressure are recorded and Fourier analyzed. Results of 
these spectra are shown in Fig. 10. 
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Mean Disturbance 

Figure 9 Contours of mean and disturbance vorticity, broadband forcing, showing probe locations 


Vorticity Specira Surface '-’ressure Spectra 




Figure 10 Spectra of disturbance vorticity and fluctuating surface pressure, at probe locations show in Fig. 9 

The broadband nature of these disturbances is again next highest on the side edge. The disturbance vorticity 

apparent; also significant is the fact that the fluctuating spectra show the expected falloff in the higher 

surface pressure is highest near the upper corner, and 
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frequencies as the disturbances convect along the shear 
layer and over the vortex. 

In the above simulations, it is assumed that the 
disturbances have no variation in the direction normal 
to the crossflow plane that is cut from the RANS 
solution. While this is a reasonable assumption for the 
instabilities that result from the cylindrical shear layer, 
it is conjectured that instabilities should exist also in the 
shear of the vortex itself, especially in light of the 
strong jet-like flow in the core of the vortex. Further, 
the disturbances resulting from this vortex instability 
should have significant oscillation in the streamwise 
direction. Assuming that the spatial scale of such 
oscillation is small compared with the scale of variation 
in the streamwise direction of the mean flow, then we 
can apply the additional ansatz that these oscillation are 
homogeneous in that direction in the simulation. For 
the following simulations, therefore, we discretize the 
three-dimensional Navier-Stokes equations, and replace 
all 3/5x terms by /a, where a is a prescribed streamwise 
wavenumber parameter. The full three-dimensional 
RANS solution, interpolated onto a crossflow plane as 
before, is taken as the mean state, and the (complex) 
equations are discretized in the crossflow plane and 
time advanced as before. 



With this assumption, the parameter space to 
be explored now has two parameters: frequency and 
streamwise wavenumber, with the results shown earlier 
corresponding to the case a=0. Results from this study 
are as yet incomplete, but it has been found for the most 
part that the shear-layer instabilities are maximally 
amplified for a=0. However, there exists a separate 
family of instabilities, associated with the conjectured 
vortex-instability mechanism; the streamwise 
wavenumber for these disturbances is fairly high, 
corresponding to wavelengths on the order of V* to Vi of 
the vortex diameter. The dominant frequency band of 
these disturbances is considerably lower than that for 
the shear-layer instability; in model frequencies, the 
shear-layer instability band was roughly 5 to 30 kHz, 
whereas the vortex instability band is about I to 10 kHz. 
In Fig. 1 1 is shown a comparison of snapshots of 
disturbance vorticity for the two disturbance modes, 
using broadband forcing for each simulation; the vortex 
mode has a streamwise wavelength of 1/3 of the vortex 
diameter. Note that the vortex disturbance has a ring- 
like structure. This is more clearly discemable in Fig. 
12, in which isolevel surfaces of disturbance vorticity 
are shown for these same modes. 



Figure 11 Disturbance vorticity snapshots, shear-layer and vortex instability modes 
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Figure 12 Isolevel surfaces of disturbance vorticity, shear-layer and vortex instability modes 


These vortex disturbance modes have 
implication in the generation of noise, for flap 
configurations for which the vortex remains in 
proximity to the flap surface all of the way to the 
trailing edge; as noted above, realistic flap sections 
frequently show this behavior. As these disturbances 
convect past the trailing edge, they would be expected 
to have considerable amplitude and be quite coherent in 
the spanwise direction; thus it would be expected that a 
significant amount of locally-enhanced trailing-edge 
noise would result. While initial experimental 
observations indicate that this is indeed the case, 
comparisons with experimental data are ongoing. 

Conclusions 

Demonstrated in this paper is a computational 
framework to predict the spectral content of fluctuations 
in the flowfield near the side edge of a flap in a multi- 
element high-lift system. The framework is 
heuristically justified by comparison with high- 
resolution steady RANS solutions of the complex multi- 
vortex system; these RANS solutions are then utilized 
as mean flowfields within which unsteady disturbances 
are simulated* Simulations show two basic families of 
disturbance modes; the first is associated with instability 
of the cylindrical shear layer, which overlies the side- 
edge separation in upstream stations, and feeds the 
trailing vortex in downstream stations. The second 
family of disturbance modes is associated with 
instabilities of the vortex and its jet-like core flow, and 
possesses significant oscillatory structure in the 
streamwise direction. 
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